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Abstract 

We derive radiative transport equations for solutions of a Schrodinger equation in a periodic 
structure with small random inhomogeneities. We use systematically the Wigner transform and 
the Bloch wave expansion. The streaming part of the radiative transport equations is determined 
entirely by the Bloch spectrum, while the scattering part by the random fluctuations. 

Keywords. Radiative transport, waves in random media, Wigner distribution, semiclassical 
limits, Bloch waves. 

1 Introduction 

Radiative transport equations describe propagation of the phase space energy density of high fre- 
quency waves in a medium with weak random impurities whose correlation length is comparable 



to the wave length A IS]. The background medium may vary but only on scales that are much 
larger than A. The phase space energy density has also been studied for a periodic potential when 
the period is comparable to the wave length [10, [D|]. It is then found that the limiting phase 



space energy density satisfies a system of decoupled Liouville equations with Hamiltonians given 
by the Bloch eigenvalues. Thus, we expect that the addition of small random fluctuations to the 
periodic structure will give rise to a system of coupled radiative transport equations. This is what 
we show in this paper formally, using asymptotic expansions. The rigorous proof of our final re- 
sult remains an open problem. Our main result is the system of the radiative transport equations 
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|) that describes the propagation of phase space energy densities a m (t, x, p) corresponding to 
various Bloch eigenvalues E m (p). We derive here the transport equations only for the Schrodinger 
equation. However, the generalization to more general types of waves in a periodic and random 
medium, like those described by hyperbolic systems considered in fpTSf ] , is straightforward. This 
results in the replacement of the eigenvalues E m {p) in (EH) by the corresponding Bloch eigenvalues 
for the hyperbolic system under consideration. The paper is organized is follows. First, we recall 
in Section 2 some properties of the Bloch eigenfunctions and also give a formal derivation of the 
Liouville transport equations in the absence of random inhomogeneities. These results were previ- 
ously derived rigorously in |l5| and |u]], and we show the connection between our formalism and 
their results in Section 2.4. Section 3 is the main part of the paper, where we derive the radiative 
transport equations in the presence of a random potential. 

2 Waves in a periodic structure 

2.1 The Schrodinger equation 

We give here a derivation of the Liouville equation for the phase space energy density like the 
one in [19 1 . A different analysis is given in |10| , [lq |. We then adapt the analysis of [|^] to the 



periodic-random case. 

There is not a lot of mathematical work on the transport limit for the Schrodinger equation with 



random potential. We cite here the work of Martin and Emch |lq| , of Spohn [20], of Dell' Antonio 



[0] and the recent extensive study of Ho, Landau and Wilkins [11|. These papers established the 
validity of the kinetic linear transport equation for a small time T > 0. The global validity of this 
limit was proved recently by Erdos and Yau in [|j] . They treat only spatially homogeneous problems 
but it is known how to extend the analysis to the spatially inhomogeneous case (slow x-dependent 
initial data and potential) ||. A really satisfactory mathematical treatment of radiative transport 
asymptotics from random wave equations is lacking at present. 

It is convenient for us to use the usual Wigner distribution and not the Wigner band-series as 
in [10, [D|. The two formulations are, however, equivalent [15]. Let e (i, x) be the solution of the 



initial value problem 



£ 2 . ,x, 



i^ + -^-V(-) = (!) 

e (O,x) = 0»( x ). 



The initial data <^(x) is uniformly bounded in L 2 (R d ): \\<ft®\\ L 2 < C. We also assume that it is 
e-oscillatory, that is, for any test function tp £ C c (R d ): 

limsup / \ip<p^Ck)\ 2 dk — ► as R goes to +00. (2) 

e^O J\V\>R/e 

Finally we assume that the family (fP e is compact at infinity: 

limsup / \4>° £ (~x)\ 2 dx. — > as R goes to +oo. (3) 

e^O J\x\>R 

A sufficient condition for (Q) is ||V^|| < C/e. The case of e-independent initial data was studied 
in ||, where the Liouville equation for the wave amplitudes (not energies) was derived. 
The potential V(z) is periodic: 

V{x + u) = V{t), 

where vector v belongs to the period lattice L: 

L = jl>i e iK eZ |> ( 4 ) 

and ei , . . . , e^ form a basis of R d with the dual basis e J defined by 

(ej -e k ) = 2ir5 jk 

and the dual lattice L* defined by (Q) with ej replaced by eP . We denote by C the basic period 
cell of L and by B the Brillouin zone: 

B = jkGR | k is closer to /z=0 than any other point \i £ L* > . 

We define the Wigner distribution by 

W e (t,x,k)= [ -^e ik ' y ^(i,x- £ y)^(t,x). (5) 

Jr<i (2ir) a 

This definition is equivalent to its symmetric version 

jRd (2ir) a 2 2 

in the sense that W £ and W e (x, k) have the same weak limit as e — ► [Sj. We also have that 

£ e (£,x) = |</> £ (£,x)| 2 = ! dkW 6 (t,x,k)= / dkW £ (i,X,k). 



The basic properties of the Wigner distributions are reviewed in detail in [[H], |l8[ . In particular 
the weak limit VF(t,x, k) of W e (t,x, k) exists in S'(R d x R d ) and under assumptions (0) and (pi) it 
captures correctly the behavior of the energy £ £ : 



lim / dx£ £ {t,x) = dxdkW(t,x.,k). 



We deduce from (|T]) and (||) the following evolution equation for W £ (t,x, k): 






+ k • V X W £ + ^A X W £ = - V e^^Vifi) [W e (x, k - /*) - W e (x, k)] . 



/xeL* 



Here V"(/x) are the periodic Fourier coefficients of V(y): 

1 



^.l-t)-T^] c dye-^-yV(y), 



so that if V(y) is smooth we have 



and 



V(y) = £ e^VM 



(6) 



(7) 



icf^-E*-")- 



ueL 



We introduce a multiple scales expansion for W^ 



W r e (t,x,k) = W (t,x,j,k) + eW r i(t,x,pk) + 



(8) 



0) 



and assume that the leading term Wo(t, x, z, k) is periodic in the fast variable z = — . As usual, we 
replace then 



v x — v x + -v, 

e 



in (^) and rewrite it as 
dW £ 

dt 



+k- 



V x + -V z 

e 



w £ + ^(v x + -V.) • (V x + -V Z )W £ 






/xeL* 



We insert the perturbation expansion (||) into ([!(]) and get at the order e : 

cw = 0, 



(10) 



(11) 



where the skew symmetric operator C is given by 

£/(z,k) = k-V z / + ^A z /-i J2 e^' z FQu)[/(z,k-,z)-/(z,k)] 



2.2 The Bloch functions 

The eigenf unctions of C are constructed as follows. Given p G R d consider the eigenvalue problem 

-^A z *(z,p)+y(z)*(z,p) = E(p)tt(z,p) (12) 

^(z + i/,p) = e iv ' u ^{z,p), for all 1/ G L 
— (z + i/,p) = e ipl/ — (z), foralli/GL. 

This problem has a complete orthonormal basis of eigenfunctions ^^(z, p) in L 2 (C): 

(y^)=J c ^% l (z,p)^(z,p) = 5 mj 6 aP . (13) 

They are called the Bloch eigenfunctions, corresponding to the real eigenvalues E m (p) of multi- 
plicity r m . Here a = 1, . . . ,r m labels eigenfunctions inside the eigenspace. The eigenvalues E m (p) 
are L*-periodic in p and have constant finite multiplicity outside a closed subset F m of p £ i? of 
measure zero. They may be arranged Ei(p) < ^(p) < ••• < ^j'(p) < ••• with Ej(p) — ► oo as 



j — ► oo, uniformly in p []2|, 13, 22]. We consider momenta p outside the set .F m . 



The problem (12) may be rewritten in terms of periodic functions <3?(z,p) = e ' p ' a \P(z, p): 

- ~ A Z <C + y(z)d>» - * P • V Z <C + ^<C = E m (p)^(z, P ). (14) 

We differentiate ( |i~4"D with respect to pj and take the scalar product of the resulting equation and 
& m to get 

^-^=^,-z(^-,^), 



which may be rewritten as 



<W = *(*mi ^T— )• (15) 



There is no summation over m in (|l5|). 

Recall the Bloch transform of a function 4>(x) G L 2 (R d ) 

C(p)= / dz0(z)^(z,p). 

It has the following properties: 

(i) 0(x) = M EE dpC(p)^(x, P ), x G i^. 

1 ' m=l o=l 

(ii) Let (/)(x),ry(x) G L 2 (R d ), then the Plancherel formula holds: 

/ dx^(x)f?(x) = i^f dpC( P )l^ 



(iii) The mapping 4> —> <ft 1S one-to-one and onto, from L 2 (R d ) — ► (B m L 2 (B). 
We deduce from these properties the orthogonality relations: 

5 & ~ x ) = m E / d p^(x, P»(y, P) 

and 

SjmSapSperiP ~ <0 = 75T / dx ^?( X ' P)*m( x > 0)- (16) 

|-D| Ji?d 

The periodic delta function 5 per in (|lq ) is understood as follows: for any function 0(p) G C°°(B) 

</>(p) = / dq(/>(q)5 per (p - q). 
Jb 

Given any vector k G R d we may decompose it uniquely as 

k = p k + Ai k (I 7 ) 

with pk G -B and /x k G L*. We then define the z-periodic functions Q^ n (z, /z, p), /1 G L*, p G -B 
by 

g^ n (z, n, P ) = j£ |p(p+/^«( z _ y , p)^( z , P ). (is) 

Then a direct computation shows that 

CQ^ n (z, n, p) = i(E m (p) - E n (p))Q^ n (z, », p) (19) 

with fj, = /x k , p = pk, so Q^fn are eigenfunctions of £. 

2.3 The Liouville equations 

Now ( |i~9| ) implies that, for any p, ker£ is spanned by the functions Q^ m , which we denote by Q^f 
(to indicate that there is no summation over m). Then (|ll]) implies that Wo(t,x., z, k) may be 
written as 

W (t,x,z,k) = Wb(t,x,z,p + /i)= J2 <^f(<»x,p)QSf(z,Ai,p), peB./xa* (20) 

with /i. = /i k , p = pk- This defines a m , which is scalar if the eigenvalue E m (p) is simple, and is 
a matrix of size r m x r m if E m (p) has multiplicity r m > 1. We call <r m the coherence matrices in 
analogy to the non-periodic case E9|. They are defined inside the Brillouin zone p G B but it is 
convenient to extend them as functions in R , L*-periodic in p. 
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Next we look at e° terms in ([T(]). We get 

dWn 

-^ + k-V x W + iV x -V z W = -CW 1 . (21) 

We now integrate both sides of (pTI) against Q"(z, fi, p) over z in C and sum over \x. To evaluate 



the left side we note that we have after summing over /i. using (g), and using (13) and a change of 
variables y — > z — y: 

E / c pQ^W,P)Q?W,p) = / CxC ^^(z,p)^( y ,p)§J( Z ,p)^(y,p) 

= 5 m j5 a a>$f3l3>- (22) 



Further, using the change of variables and equations above and also (|15|), we get 

/J,£L* 

= E / ^^e i(p+MHyi - y2) fe + /^0^(z-yi,p)^(z,p)^(z-y 2 , P )^(z, P ) 
i Jcxc \C\- i dyi ■> J 

= ^mjSpP'i-Q^^j) = Smj^aa'Spp'-Q-T 1 - (23) 



A similar calculation shows that the third term on the left vanishes: 
{ 9Q^' ^ f dzdyd^(z-y,p)- ta 



y ,*££_ a**) = f dzdyd* m {z-y,p) ~ 

f? T } 8zi '^ ' Jc>,c \C\ 2 dzi j( y,pjw jl x, Pj w ml z, Pj 

f dzdy Ta i, . = , a, <9\M'(z,p) 



= —-O a /3Of30'Ojrn-K r -O a 0O00'Oj m — = U. 

The right side of ( |2~I| ) integrated against Q? vanishes since £ is skew symmetric, and Q" 6 ker£. 



Putting together (|2l|), (|22| ) and (23) we obtain the Liouville equations for the coherence matrices 

-^ + V p E m ■ V x a m = 0. (24) 

The initial data for equations ( p4[ ) is constructed as follows. Let W®(x., k) be the Wigner transform 
of the initial data </>g(x). Then a m (0, x, p) is given by 

with the limit understood in the weak sense. 



2.4 The Wigner band series 

The Liouville equations ( p4| ) were previously derived in [15, [jl| in terms of the Wigner series defined 
by 

We (t,x,k) = 7 H I E e4k '^( x -^)^( x )= E w £ ( x ,k + Ai ). 

It was shown that the weak limit w(t, x, k) of w £ (t,x, k) has the form 

w(£,x,k) =E^j( i > x : k )- 

Here Wj is the limit Wigner series of the projection <j>> e of the solution <p e of the Schrodinger equation 
on the Bloch spaces <SJ: 

S i = {/ e W) = /(x) = E ||| / fl J^/«(p)*?(f ,P)} • 



Each projection Wj evolves according to the Liouville equation (|24|). This result may be related to 
our approach as follows. The Wigner distribution W e (t, x, k) has the asymptotics 



H/ £ (t,x,p + A *)« E < /3 (^x,p)Q^( 7 ,/i,p). (25) 

m,a,(3 



Then the weak limit of w e is given by 

w(t,X,p)= E O'wfC^^P) E / TT^Qm ( Z ,/^P) =E Tr °' m ( f ' X 'P)- 

If we take the trace of ( pl| ) we get the transport equations for w m (t, x, p) = Tr<7 m (£, x, p) obtained 
in p5[ and [In]]. However, the representation ( |25| ) captures not only the weak limit of w e but 
also the oscillations of W s on the fine scale. Moreover, the cross-polarization of various modes 
corresponding to the same eigenvalue is also taken into account by the off-diagonal terms in the 
coherence matrices. The oscillations of the spatial energy density, which is the physically interesting 
quantity, are given by 

£ e (t,x,p)«E / d P £7 m(^x,p)*°(-,p)f^(-,p)^E / dpTVa m (t,x,p) ase^O. 
m,a,p m 

This information on the energy oscillations may be useful in numerical simulations. We see that 
a^ 1 are phase space resolved energy densities of different modes inside the Brillouin zone. 



3 The random perturbation 

3.1 Simple eigenvalues 

We assume first that the Bloch eigenvalues E m (p) have multiplicity one for all p G B. This 
assumption is known to be true for the leading eigenvalue when the Fourier transform (pi) of the 
periodic potential V(y) is negative Q. In many physical problems the absence of level crossings 
restricts our results to the lower part of the spectrum. We consider now small random perturbations 
of the periodic problem (|l]) with randomness being on the same scale as the periodic potential but 
weak: 

Here N{y) is a time independent mean zero spatially homogeneous random process with covariance 
tensor R(x) defined by: 

(N(y)N(y + x)> = i?(x), (iV(p)iV(q)) = (27r) d R( q )8(p + q). (26) 

Here (•) denotes the ensemble average and the Fourier transform iV(q) is 

iV(q) = / dxe- iq * x iV(x). (27) 

JR d 

The Wigner distribution W e (t,x, k) satisfies the evolution equation 

^ + k-V x W £ + ^A x W £ = rE^^W^-/*)-^)] 
m l l£ /xeL* 



V e JR d K 27r ) 



Here V"(/Lt) is the periodic Fourier transform (0) and -/V(q) is the Fourier transform (|27| ) over i? rf . 
We consider the asymptotic expansion 

W e (t,x,k) = W (t,x, *,k) + >/eWi(t,x, *,k) + £W 2 (t,x, *,k) + ..., 

with the leading order term Wo being deterministic. We introduce as before the fast variable z = — , 

£ 

replace V x — ► V x -\ — V z and collect the powers of e. The order e _1 gives as before 

e 

£Wn = 0. 



Thus we still have the decomposition (|20|): 

W {t, x, z, p + n) = Y^ Vm(t, x, p)Q m (z, n, p) (28) 

m 

with a m being scalar because the spectrum is simple. Recall that the functions a m (t, x, p) are 
L*-periodic in p and the z-periodic functions Q m are given by (^) with m = n. The order e~ 1 ' 2 
terms give 

£Wl = ~ I 7% ei * Z ti(<l) [Wb(z, k - q) - W (z, k)] . (29) 

The distribution W\ need not be periodic in the fast variable z. Therefore it may not be expanded 
in Q mn and we use the the basis functions 

P mn (z, v, P, q) = J c 0re^P + ^> y * m ( Z - y, p)fc m (z, p + q), (30) 

defined for z G i? rf and p, q£B, in place of the periodic functions Q mn (z, /z,p). The functions 
P mn are quasi-periodic in z: 

-Pmn(z + V, /i, P, q) = -Pmn(z, /X, p, q)e~ lUfi . 



We use ( |l6|) and (||) to obtain the orthogonality relation for the functions P, 

r dz 
IriJb 



r dz 

Yl / TBi p mn{2,fJ-,-p,q)Pji{z,fi,p,q ) =S mj 5 n i5 P er(q.-qo)- (31) 



The operator £ acts on these functions as 

C-Pmn = i(E m (p) ~ E n (p + q))P mn (z, /x, p, q), 



which is the analog of ( fL9[ ) for the functions Qmn- Note that the integrand in ( p0[ ) is periodic in y 
so no boundary terms are produced by integration by parts. 
We decompose W\ in this basis as 

Wi(t,x,z,p + n) = J2 r^j-^mn(i,x,p,q)P m „,(z,/i,p,q), (32) 

with z G P d , p£B and /x G L*. We insert (p 7 ^) into (p9[), multiply ( |29| ) by P/^z, /x,p,qo), sum 
over //£l* and integrate over z G P . Then we get using ( |3l|) on the left side 

,, \ ST ff dzd<ie l * z N(ci) [W (z, p + /a - q) - V^ (z, p + /x)] - 
r ?i Kt,x,p,q )= L /L, T^Ad F .^^_ F .^^fl ^(z,/x,p,q ), (33) 



/lG /» 7 ^ (2n) d Ei(p + q ) - ^ (p) + *0 
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where 9 is the regularization parameter. We let 6 — > at the end. We insert expression ( |28| ) for 
Wo and the definition ( |TS| ) of Q m into fl33|). The resulting expression may be simplified using (O) 
and (§): 

^(t,x, P , qo ) = (^ E El{p + qo) _ Ej{p) + ig ^-(p + qo + /*, P) 

rfzrf q c »q. z ^(q)^- (p) ^ ( z > p + qo)^j ( z > p) / 34 x 

fld XjRd (27r) d |B| Ei(p + q )-E j {p)+id l j 



Here the amplitude A\j (q, p) is given by 



4i(q,P) = / c ^e-^- p >^Ky,q)^(y,p). (35) 



The next order equation is 

dW . _ TTr _ _ TTr „ TT7 I /" dq 



, i; +k-V x Wo + tV z -V x Wo + £W 2 = - I ? -%e it i- z iV(q)[Wi( Z ,k-q)-Wi(z ) k)]. 
We multiply this equation by Qj(z, /x, p), integrate over z and sum over fj, £ L*, and take average. 



Then as before the left side is 



doj_ 

at 



LHS = -^- + V p Ej • V x CTj. (36) 



The right side is 



RHS = /i + / 2 , (37) 



where 



^iEj c p/ B ^W«P + ^^AM)) (38) 



/xeL* 



and 



We insert expression ( |32| ) for V^i into (|3^) to get 

h = - J2 / tftt / 7^~w / T7iT^( c i) ( 5i( z '/ x 'P)I] 7 ?™™(p + / i ~q'qo)- p ™™( z '/ i 'P~q'qo))- 

We may split 7i = I\\ — I\% according to the two terms in (^). We insert the first term in (p^) 
into the expression for In and average using spatial homogeneity ( P6|) of the random process N(z), 
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orthogonality (|l3|) of the Bloch functions x f m (z, p), and also sum over /x £ L* using (||). Then we 

get 

^ 1 y> ^ /• rfq i?(q + /i)|^ im (p,p-q-/j)| 2 a j (p) 
11 i ^l^^\B\{2^Y Ej(p)-E m (p-^+ie 

Here we have replaced integration over q £ i? d by integration over £> and sum over /x £ L*. The 
second term I12 is evaluated similarly: 

_ 1 f dq J?(q + /i)|^ im (p,p-q-/x)| 2 o- m (p-q) 

12 i ^L*mJB\B\{2-K) d E:i (p)-E m (p- q ) + i9 

Thus we have 

r I V V /" dq -R(q + /^)l^im(p,p-q-M)| 2 [q-j(p) -q- m (p-g)] , , 

1 V^ -^ |B|(27r) " ^■(p)- J B m (p-q)+^ • l j 

One may verify that I2 = ^1 in (|37|)- We insert then (|39| ) into (|37|), take the limit # ^ 0, make a 
change of variables q — > p — q, and combine it with (36) to get the system of radiative transport 
equations: 

§ + V p Ej • VxCT, = J2 I T§iQjm(P, q)ktn(q) - ^(pM^-Cp) - S ro (q)). (40) 

at - Jb \B\ 

The differential scattering cross-sections <2j m (p,q) are given by 

Gjm (p,q)= I] /„ w _! ^(p-q + M)l^jm(p + M,q)| 2 - 

This is the main result of this paper: we have derived a system of coupled radiative transport 
equations for the phase space energy densities of the Bloch modes. The transition probabilities 
Qjm(p>q) are real and symmetric: Qj m (p,c{) = Qmji^P) as seen from the definition ( p5[) of 
Aij(p, q). Therefore the total energy is conserved: 



£(t) = 22 / dx. I dpa m (t, x,p) = const. 



Transport equations like (^) are well known in the theory of resistance of metals and alloys [17|. 
Their systematic derivation from the Schrodinger equation with a periodic and random potential 
(|26l) is new. 



A Appendix. Multiple eigenvalues 

When eigenvalues Ej(p) are not simple, but their multiplicity is independent of p, so that there 
are still no level crossings, the analysis in the previous section can be extended to this case. This 
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case is probably very rare for the Schrodinger equation but is important for other types of waves, 
for instance, in symmetric hyperbolic systems. We present it here for the sake of completeness. 
The result is as follows. Let the matrices T°^(p, q), a = l,...,tj, (3 = 1, . . . , r m , where Tj and r m 
are the multiplicities of Ej and E m , be defined by 

2#(P,q)= / f e^-q)- z ^(z,q)^(z,p). 

JC7(27T) — \C\ 

Then the coherence matrices CTj(p) satisfy the system of radiative transport equations 
dan 

■ + v p 

dq 



-£ + V v E r V*aj 



£ I ^i?(p-q + /x)r jm (p,q-/i)a m (q)T* n (p,q-/x)5(^(p)- J E m (q)) 
TT* JB -D 



m£JV 



i f <iqi?(p — q + fx) 



2vr J B IB] 



Tjm(p, q - v)T* m (p, q - /i)oj(p) o- j (p)T i?n (p, q - /x)T/ m (p, q - /x) 



£? J -(p)- J E7 m (q)+iO %(p)-.E m (q)-iO 

These equations have the same structure as the radiative transport equations for polarized waves 
derived in |01]. The expression of the total scattering cross-section as a principal value integral 
and not as a familiar integral against 5[E m (c[) — Ej(p)] is known in transport theory for polarized 
waves |l2| , Q PL9f1 . They reduce typically to the form common in scalar transport equations under 
additional symmetries, like rotational invariance of the original wave equations and the power 
spectrum tensor. 
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